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Abstract. 

Properties of the four families of recently introduced special functions of two real variables, denoted 
here by E ± , and cos ± , are studied. The superscripts + and ~ refer to the symmetric and antisymmetric 
functions respectively. The functions are considered in all details required for their exploitation in Fourier 
expansions of digital data, sampled on square grids of any density and for general position of the grid 
in the real plane relative to the lattice defined by the underlying group theory. Quality of continuous 
interpolation, resulting from the discrete expansions, is studied, exemplified and compared for some 
model functions. 
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1. Introduction 

A large body of empirical evidence as well as theoretical experience has been accumulated in treating 
two-dimensional digital data owing to the large amount of data requiring treatment in practical contexts. 
We approach the problem from the opposite direction. Our departure point, rather than a set of specific 
2D data, are properties of new special functions in n dimensions [6l[7]. This presents certain advantages 
and disadvantages, some of which are listed in the concluding remarks. The most important distinction 
is undoubtedly the possibility to carry out any analysis of lattice data in the Fourier space rather than 
in the data space. 

In this paper, common cosine transforms in one and two dimensions [T8| IT3] are based on the sym- 
metric group S2 and S2 x S2 respectively. Our approach to 2D problems is built on group S3 and on 
its afflne extension. The special functions differ essentially from the product of trigonometric functions 
each depending on a single variable extending in mutually orthogonal directions in the case of S2 x S2. 
We start with the properties of the special functions of S n [6, 7J. 

Our aim here is to take advantage of the fact that, when working in 2D, many specific properties of 
functions can be expressed explicitly and more detailed questions can be answered than in the case of 
the general dimension [HI [7] . More precisely, we consider two families of symmetric and antisymmetric 
special functions of two real variables. They can be viewed as generalizations of the common exponential 
and cosine functions of one variable. 

For each family, we study: 

(i) The orthogonality of the functions on a lattice of chosen density within an appropriate finite region 
of the real Euclidean space M 2 , and the corresponding finite Fourier expansions of digital data functions 
sampled on such lattices. 

(ii) The interpolation technique of the digital data used for the expansions into a continuous differ- 
entiable function. Interpolation quality is assessed as a function of the density of the lattice. Cosine 
transforms of types I, II, III, and IV [18] can be generalized by situating the points of the sampling grid 
in a general position within the orthogonality region. It is important to point out those versions of the 
formalism by which the undesirable Gibbs effect can be avoided |17j . 

The definition of the symmetric and antisymmetric exponential functions E + and E~ is analogous 
to the definition of so called C— and S— functions, where instead of the symmetric group 5 n , the Weyl 
group W of some compact semisimple Lie groups is used [TQlUT]. So called E— functions [12] which are 
based on even subgroup of the Weyl group W are analogous to 'alternating' exponential functions [3 [9]. 
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In some cases, a direct relation between these two concepts can be found [16J. The coincidences of the 
orbit functions are due to the known isomorphism of the symmetric group S n and the Weyl group of 
the Lie group SU(n), and of the alternating subgroup of S n and the even subgroup of the Weyl group. 

In Section 2, the antisymmetric and symmetric 2-dimensional functions are defined respectively as 
2x2 determinants and permanents [15] of exponential functions of one variable. Properties of the 
antisymmetric functions (determinants) are considered first, namely their continuous and discrete or- 
thogonality, the discrete Fourier transforms, and others. An exposition of analogous properties of the 
symmetric functions (permanents) follows. Examples of real and imaginary parts of the functions are 
shown in Fig. 1 and 2. 

Section 3 is devoted to 2D interpolation by symmetric and antisymmetric exponential functions. 
General 2D interpolations are recalled, followed by a detailed study of interpolation by antisymmetric 
functions. Examples of the interpolation of the model function in Fig. 3 are shown in Fig. 4 for several 
densities of the sampling grid. Interpolation properties of the symmetric exponential function follow. 
Examples of the interpolation of the same model function are shown in Fig. 5. Interpolation errors are 
summarized in Table 1. 

Section 4 parallels Sec, 2 and 3, except that the exponential functions are replaced by 2 x 2 antisym- 
metric (determinants) and symmetric (permanents) of cosine functions. Examples of antisymmetric and 
symmetric cosine functions are shown respectively in Fig. 6 and 8; interpolation examples can be found 
in Fig. 7 and 9. Interpolation quality in the above examples is compared in Table 1. The difference 
between the model function and its interpolations is integrated over the whole orthogonality region. 

The last example deals with known undesirable effect of Fourier interpolation, referred to as the 
Gibbs effect [17J. In particular, it is shown that the Gibbs effect is absent in two versions of our 2D 
cosine transforms, as illustrated in Fig. 10. 

The last section contains various related comments and remarks. 

2. Two-dimensional ( antisymmetric exponential functions 
2.1. Two— dimensional antisymmetric exponential functions. 

2.1.1. Definitions, symmetries and general properties. 

Two-dimensional antisymmetric exponential functions E^ x ^ : R 2 — > C have for A, fi G K the following 
explicit form 



p 27ri\x p 27ri\y 



27ri(\x+ny) _ 27ri(\y+nx) 



jp- (rr n,\ — e^ %AX e^ tA y _ p 27ri{Ax-\-fiy) _ 2m{Ay+tj,x) /i\ 



We observe that E^ xx ^(x,y) = and E^ x ^{x,x) = 0. From the explicit formula we immediately 
obtain antisymmetry of E~^^(x,y) with respect to the permutation of variables y) 

E (x,^(y^) = - E (x,^y)- ( 2 ) 

and, moreover, with respect to permutation of (A, fi) 

E (x,^y) = -%,x)(^y)- ( 3 ) 

Therefore, we consider only such E^ x ^ with so called strictly dominant (A,/x), that is, pairs (A, /x) 
with A > /i. 

The functions E^ k ^ with fc, Z E Z have additional symmetries related to the periodicity of exponential 
function. One can verify directly that 

E (k,i)( x + r >y + s ) = E (k,l)( x ^y^ r,s eZ (4) 
The relations Q and (jij) imply that it is sufficient to consider the functions E~^ kl y k,l E Z on 

the so called fundamental domain F(S r | ff ) [6j. The following open 'triangle' can be chosen as the 
fundamental domain in dimension two, 

F(Sf) = {(x, y) G (0, 1) x (0, 1) | x > y} . (5) 
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We also have the following property for a E R 

E^(x + a,y + a) = e 2 ^ k ^ a E^ l} (x,y). (6) 

2.1.2. Continuous orthogonality. 

The functions E~ are mutually orthogonal on the fundamental domain, i.e. 

/ E~Jx,y)E7 yi Jx,y)dxdy = 5 k k'$ii', k,l,k',l' E Z, k > l,k' > l'. 

JF(S<f) K ' ; V ' ; 

Therefore, every function / : R 2 — > C that is antisymmetric /(x, y) = — /(y, x) and periodic f(x + r, y + 
5) = /(x, y), r, 5 E Z and has continuous derivatives can be expanded in the antisymmetric exponential 
functions E7 



J (k,iy 



f(x,y)= V c k iE (kl) (x,y), c M = / f(x,y)E (kl) (x,y) dxdy. 



I fe>Z ; 

2.1.3. Solutions of the Laplace equation. 

The functions are solutions of the Laplace equation 



(7) 



{& + £2) E (k,l)(^y) = -^ 2 (k 2 + l 2 )E- kl) (x,y) 



(8) 



and, moreover, of the equation 

d 2 d 2 



j^^2 E (k,l)( x >y) = l^ 4 k 2 l 2 E (kl) (x,y) 



which is algebraically independent. The functions E satisfy the condition E, kl Jx,y) = on the 
boundary x — y. 

2.1 A. Discrete orthogonality. 

Discrete orthogonality of antisymmetric exponential functions over the grid of form (x m ,y n ) = 
(m/N, n/N), m,n E {0, . . . , N — 1}, m > n was proved in [6J. The positive integer N sets the density 
of the grid inside F(S^). For applications, it may be convenient to consider orthogonality over more 
a general type of grid. In addition to the parameter TV, we parameterize the grid by parameters a E R 
and b E [0, 1]. The equidistant N(N — l)/2-point grid L~ bN1 is given by 

L a,b,N,l i( X mi Dn) I ™, > 71, 777, U = . . . N - 1} 

where 

/ \ f m + b n J r b\ 

(x m , y n ) = ( a + + I . (9) 

Using the property Q, we observe that the orthogonality relations from [6J also hold over the grid 

L a,b,N,V 

N-l 



^2 E (k,i)( x ™> yn)E( k ,j,)(x m , y n ) = N 2 5 kk f5 w , (10) 

f m,n=0 1 
\ m>n J 

where fc, /, A; 7 , Z 7 E {0, . . . , TV - 1}, k> l,k f > V . The set of functions ^ ^ with fc, Z E {0, . . . , 3}, k > I of 
mutually orthogonal functions on the grid L~ hA 1 is depicted in Figure [Li 

Note that for N odd, iV = 2M + 1, the discrete orthogonality relations (10) also hold for the set of 
functions E^ k ^ with fc, / E {-M, . . . , M}, fc > Z. 



JIRI HRIVNAK 1 ' 2 AND JIRJ PATERA 1 






real part 



imaginary part 






real part 



imaginary part 



Figure 1. The contour plots of continuous functions n, k,l E {0, ...,3},fc > /, 
which are discretely pairwise orthogonal on the grid L~ bA r 



2.1.5. Antisymmetric discrete Fourier transform. 

Suppose we have a discrete function / : L~ b N 1 — > C defined on the grid L~ b NV The antisymmetric 

discrete Fourier transform of / over L~ b N1 is given by 



JV-l 



@ki = J[2 52 f( x ™,yn)E {k l) (x m ,y n ) k > I, fc, I = . . . N - 1. 



(11) 



f m,n=0 1 
\ m>n J 

Orthogonality relation (|10|) immediately gives the inverse transform of N(N — l)/2 coefficients /?^: 



AT-l 



ffc,Z=0\ 
I fc>Z J 



(12) 



2.2. Two— dimensional symmetric exponential functions. 

Two-dimensional symmetric exponential functions E^ x , : R 2 — > C have for A,/iGK the following 
explicit form 



(13) 



From the explicit formula we immediately obtain the symmetry of E~^^(x,y) with respect to permu- 
tation of variables (x, y) 
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and, moreover, with respect to the permutation of (A, ji) 

Therefore, we consider only such Ef x ^ with so called dominant (A, /x), that is, pairs (A, jn) with A > \i. 
Functions with k,l E Z have additional symmetries related to the periodicity of exponential 

function. One can verify directly that 

E (kj)( xJrr 'y + 5 ) = E (k,i)( x 'y)i r,s (16) 



The relations (|14|) and (|16|) imply that it is sufficient to consider the functions E^iy M ^ Z on the 



closure of the fundamental domain F(5| ff ) [6]. 
We also have the following property for a E I 



E+ kl) (x + a,y + a) = e 2 ^ k+l ^ E+^x , y) . (17) 



2.2.1. Continuous orthogonality. 



The functions E + are mutually orthogonal on Ff^S^), i.e. 



J ^ ^E^ kl) {x,y)E^ kll ^{x,y)dxdy = Guhk'^lV , \l' E Z, k>l,k f >l' 



where symbol G k i is defined by 



Gu= 2 (18) 
1 otherwise. 



Every function / : R 2 — > C that is symmetric f(x,y) = f(y,x) and periodic /(x + r, y + s) = 
f(x,y), r, 5 E Z and has continuous derivatives can be expanded in the antisymmetric exponential 
functions E^~ k ^ : 

f(x,y)= V c kl Et kl) (x,y), c kl = G~} \ f(x,y)E+ kl) (x,y)dxdy. (19) 

2.2.2. Solutions of the Laplace equation. 

The functions are solutions of the Laplace equation 

(£> + ^) 4,0^) = -4vr 2 (fe 2 + I 2 ) J E+ (20) 



and moreover of the equation 



^2^ E tk,l)^y) = 16vr 4 fc 2 ^ 2 ^+ (x,y). 
The functions E + satisfy the condition 

where n is the normal to the boundary x — y. 

2.2.3. Discrete orthogonality. 

Discrete orthogonality of symmetric exponential functions over the grid of form (x m , y n ) — (m/N,n/N), 
m,n E {0, . . . , TV — 1}, m > n was proved in jB]. We consider the orthogonality over a more general 
type of grid. Besides the parameter TV, we parameterize the grid by parameters a E R and b E [0, 1]. 
The equidistant N(N + l)/2-point grid L^ b N1 is given by 

L tb,N,i {( x ™, y n )\m>n,m,n = 0...N-l} 
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where (x m ,y n ) are given by Q. Using the property (17), we observe that the orthogonality relations 
from [6J also hold over the grid L^ bN1 : 



N-l 



^2 G mn^,Z)(^m,yn)^/^)(^m ? 2/ri) = G ki N 2 5 kk > 5 W , 
f m,n=0 1 
I m>n J 



(21) 



where fc, /, A;', /' G {0, . . . , TV - 1}, k>l,k f > V \ The set of functions £7+ ^ with fc, / G {0, 1, 2}, fc > / of 
mutually orthogonal functions on the grid L~^ b 3 1 is depicted in Figure jlj 

Note that for TV odd, TV = 2M + 1, the discrete orthogonality relations ( |2l| ) also hold for the set of 
functions Et z) with fc, / G {-M, . . . , M}, fc > /. 






real part 



imaginary part 






real part 



imaginary part 



Figure 2. The contour plots of the continuous functions E^ kl y k,l E {0,1, 2}, A: > /, 
which are discretely pairwise orthogonal on the grid L^ b3 r 



2.2.4. Symmetric discrete Fourier transform. 

Suppose we have a discrete function / : L^ bN1 —> C defined on the grid L^ bN1 . The symmetric 

discrete Fourier transform of / over L^ b N 1 is given by 



N-l 



Pkl=Q^2 E fc>/, fc,/ = 0...TV-l. (22) 

kl f ra,n=0 1 
1 m>n J 
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The orthogonality relation (10) immediately gives the inverse transform of N(N + l)/2 coefficients (3^: 

N-l 

f(Xm,yn)= ^ Pkl E (k,l)( X ™>yn) ( 23 ) 
(k,l=0\ 
I k>l J 

3. Interpolation by ( antisymmetric exponential functions 

3.1. General two— dimensional trigonometric interpolation. 

Let us consider a symmetrically placed square in R 2 with a side of T G 1. For aEi, the square is 
given by 

K[ a ,a>} = [a, a] x [a, a] (24) 
Let us also select an arbitrary natural number N and parameter b E [0, 1] and consider a symmetrically 
placed 7V 2 -point grid L«,&,tv,t = {Om, Vn) \ m, n = . . . N - 1} C K [a ^ T] where 

Om, 2/n) = ( a + ~^^ T ' a + ^\^ T 

Suppose we have a given function / : if[ a?a +T] ~~ ^ C and a set of points L a ,6,A^,T C ^[ a ,a+T]- I n the 
following, we distinguish two cases, namely N — 2M+1, or TV = 2M. The (trigonometric) interpolation 
problem can be formulated in the following way: find a trigonometric interpolating polynomial of 
the form 

M 

t/>N 9 Tfav) = Ckie 27rik ^e 2 * il T, x,y eR (25) 

k,l=-M 

such that it coincides with / on the grid £ a ,6,7V,T 5 which means satisfying for all (x m ,y n ) E L a ^^,T 
the condition ^N,T(x m ,y n ) = f(x m ,y n ). Henceforward we set for simplicity T = 1, i.e. we have the 
trigonometric interpolating polynomial = ^iv,i °f the form 

M 

t/> N fav) = c kl e^ ikx e 2 ^y : x,y eR (26) 

k,l=-M 

satisfying on K [a ^ a+1] 

^N{x m , y n ) f(xm, Vn), m, n = . . . N - 1. (27) 
Note that, in all of the following formulas, we can always recover an arbitrary size T simply by linear 
transformation 

fay) (^^) • 

For N = 2M + 1, the trigonometric interpolating polynomial has (2M + l) 2 = TV 2 unknown 
coefficients c^, to which correspond TV 2 constraints (27). For N = 2M, denoting tn = e 27rl ( Na + b \ we 
assume further 4M + 1 conditions 

Ck-M = / r N Ck,Mi k = -M ...M -1 

C- M ,i = t n c m ,u l = ~M...M (28) 

and we have (2M + l) 2 - (4M + 1) = (2M) 2 = TV 2 unknown coefficients cm corresponding to TV 2 
constraints ([27|). 

aalyze tl 

/ e 27ri(-M)x e 27ri(-M+l)x ... i ... e 27ri(M-l)x e 27Ti(M)x \ 

27ri(-M)m e 27ri(-M+l)xi ... j ... e 27ri(M-l)m e 27ri(M)x : 



In order to analyze the system of linear equations ( |27| ), we introduce an N x (2M + 1) matrix 
Vm(xo, • • • ,xn-i) ' 



y e 27Ti(-M)x N _ 1 e 27ri(-M+l)xjv-i ... 1 ... e 27ri(M-l)xjv-i e 27ri(M)xjv-i y 



and the product 



N-l 

27T\X n 



W N = W N (x , xn-i) = W N (y , y N -i) ■= Yl (e 2iriXm ~ e 

f m,n=0 1 
\ m>n J 
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The coefficient matrix of the system ( |27[ ) can be written for N = 2M + 1 as 

Vm(xq, • • • , x N -i) ® Vm(vo, • • • , Vn-i)- (29) 



The matrix (29) is Vandermonde-like and its determinant is 



det[y M (^o, • • • , x N -i) ® V M (yo, • • • , Vn-i)] = t n N{N 1} W% N . 



For TV = 2M, we add to the matrix (|29j) an additional block of rows (28). The determinant of the 
resulting matrix is equal to 

aN -N 2 +2N+1 W 2N 
-4 T N W N . 

Note that both determinants are always non-zero. Thus, the solution of the interpolation problem 
always exists and is unique. 

The coefficients c^i are given for N — 2M + 1 by 

N-l 

N 2 

m,n=0 

and, introducing the symbol gu,M by 

\\ if k = -M,M 

9k,M = y, . (3i) 

I 1 otherwise, 

one can write for N — 2M 

N-l 



<7, - ' f(^,y n )e- 2nikx -e- 2 ^ (30) 



Ckl 



9k,M9l,M 



N 2 

m,n=0 



E f^m,y n )e- 27rikx -e- 2 ^. (32) 



3.2. Antisymmetric interpolation. 

For interpolation with antisymmetric exponential functions, we consider the triangle a+1 j inside 
the square K^ a ^ : 

K {a,a+i] = ?/) e [a, a + 1] x [a, a + 1] | z > y} . (33) 

For a given function / : a+1 j — > C and a set of points N 1 C a+1 j we define an antisymmetric 
trigonometric interpolating function 

M 

^-(x,y)= c ki E (k,i)( x ^^ x.yeR (34) 

( kJ=-M\ 
I fc>Z J 

satisfying 

^(%, J/n) 2/n), ra > n, ra, n = . . . TV - 1. (35) 

For TV = 2M + 1, the antisymmetric interpolating function has N(N — l)/2 unknown coefficients 
c/e/, to which correspond TV(TV — l)/2 constraints (35). For TV = 2M, we assume further 2M conditions 

cf_ M ~ t n c m,v 1 ~ M + 1 • • • M ~ 1 



C M-M — 



(36) 



and we have (2M + l)2M/2 - 2M = TV(TV - 1) /2 unknown coefficients corresponding to TV(TV - 1) /2 
constraints ([35]). 



Proposition 3.1. There exists a unique antisymmetric interpolating function (34) satisfying (35). The 
coefficients are given for TV — 2M + 1 by 

I N ~ l 

C kl ~ ~jy2 f( x miyn)E(hi^( x rn,yn) (37) 

f m,n=0 \ 
\ m>n J 
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and for N = 2M assuming ( 36 ) by 



9k,M9l,M 



N-l 



-kl 



f\j2 f( x miyn)E^ k ^(x m ,y n ). 

f m,n=0 1 
\ m>n J 



(38) 



Proof. For the function / : K^ a a+1 j — » C, we define its antisymmetric extension Af : if[ a?a +i] — > C by 
the formula 

' f(x,y) x>y 
Af(x,y) = lo x = y (39) 

K -f(y,x) x<y. 

We obtain a unique trigonometric interpolating polynomial Aifiw of the form 

M 

A^ N (x,y) = J2 (Ac)kie 2nikx e 2 ^y, x,yeR (40) 

k,l=-M 

satisfying 

A^ N {x m , y n ) Af(x m , y n ), m, n = . . . N - 1. (41) 
According to (30), the coefficients (Ac)^ are given for TV = 2M + 1 by 

N-l 



(Ac)kl = ^ £ A/(x m ,y n )e- 27ri ^e- 27ri ^. 



(42) 



m,n=0 

Note that since Af is antisymmetric we have (Ac)ki = —(Ac)ik, and so, (Ac)kk — 0. Taking this into 
account, the formula (40) can be rewritten as 

M N-l 
A^ N (x, y)= J2 (Ac)kie 2nikx e 2nily - (Ac)i k e 27rikx e 27rily 

I fc>Z J I k<l J 

M ^ 



r k,i=-M\ 



Using the symmetricity of the interpolating grid x m = y m and the antisymmetricity of Af we rewrite 



(42) 



N-l 



N-l 



(Ac) kl =^ E A f( x ^ y n )e- 2 ^ kx -e- 2 ^ - -L E /(f»> ^K^^e" 2 ^" 



f m,n=0 1 
\ m>n J 



AT2 



f m,n=0 1 
\ m<n J 



JV-1 



(44) 



~~ TV 2 ^ f( x m,yn)E^ kl ^(x m ,y n ) :— c fc/ 

f m,n=0 1 
I m>n J 

and analogously for = 2M. 

The existence of two different antisymmetric interpolation polynomials would imply the existence of 
two different interpolating polynomials satisfying (40), (41) — a contradiction. □ 

3.2.1. Calculation of the coefficients c^. 

Instead of the direct calculation of the coefficients c^, the antisymmetric discrete Fourier transform 
11) can be used, and the resulting coefficients (3^ transformed to c^'s. By direct comparison of (37), 
38) to (fTTl), we obtain for N = 2M + 1 



c kl = (3 w fc>Z,M = 0...M 
\-i = -tnPm-i^ k = 0...M,l = l...M 

2 



(45) 



-k-l ~ T N^N-k,N-V k <l, k,l 



1...M 



10 
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and for N = 2M 



c kl gk,M9i,M@kl> k > I, k,l = . . . M 
\-l ~ ~9k,M T N^N-l^ = . . . M, 1 = 1. . .M -1 
C -k-l = T N^N-k,N-V k <l, k,l = l...M -1 



(46) 



The formula (46) determines (2M)(2M — 1)/2 coefficients; the rest of the 2M coefficients are determined 
via the relations (36). 



3.2.2. Trigonometric form of i(j^(x,y). 

The antisymmetric interpolating polynomial (34) can be brought to its 'trigonometric form'. Intro- 
ducing the symbol h k via the relation 



f i if k = 
hk = S 

1 otherwise 



(47) 



we have 



where 



ip N (x,y) = hkhi[A kl (cos27vkx cos27vly — cos 2tv I x cos2ty ky) 

fk,i=o\ 

I fe>Z J 

+5^ (sin 27T fcx cos 2nly — cos 27r/x sin 2nky) + 

(cos 27rfcx sin 2nly — sin 27r/x cos 2nky) 
+D kl (sin 27rfcx sin 2nly — sin 27r/x sin 2nky)] 



kl k 5 Z ^ /c 5 ^ ^ fc 5 ^ ^ fc 5 ^ 
-^fcZ =i ( c fc,Z _ c -k,l + c fc,-i _ c -k,-l) 

c h =i (Cfc,/ + c_fc,/ - - c: fc _ z ) 

Substituting equations (37), (38) into (49) we obtain the following explicit formulas for TV 
f 2\ 2 N ~ 1 

A kl = I — J ^2 f( x ™, Vn)(cos 27vkx m cos 2nly n - cos 2nlx m cos 27vky n ) 

\ m>n J 



(48) 



(49) 



2M + 1 



B 



ki 



2 N-l 

M I ^2 /(^m, ?/n)(sin 27rfcx m cos 2nly n — cos 27vlx m sin 27vky n 

f m,n=0 1 
\ m>n J 

2 JV-1 

/V / X] f(x rn ,y n )(cos27Tkx rn sm27Tly n - sin 27r/x m cos 27r% n 

f m,n=0 1 
I m>n J 

/2\ 2 N_1 

^ V7v) ^ /(x m ,^ n )(sin27r/cx m sin27r/?/ n - sin27r/x m sin27r% n 

^ ' f m,n=0 1 
\ m>n J 



c 



(50) 
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and for N = 2M 



2 JV-1 



A kl =gk,M9i,M [jyj ^2 f(x m , y n )(cos 27vkx m cos 2nly n - cos 27vlx m cos 2nky n 



f m,n=0 1 
\ m>n J 

2 iV-1 



=gk,M9l,M ( ^ ) /(^m,yn)(sin27rA:x m cos27r/y ri - cos 2yr/x m sin 2tt/c^ 



Cm =gk,M9i,M 



f m,n=0 1 
\ m>n J 

2 JV-1 



/v r / 5^ /(x m , y n )(cos 2irkx m sin 27rZy n — sin 27ilx rn cos 2iihy n ) 

f m,n=0 1 
\ m>n J 



2 JV-1 



£^ =9k,M9i,M ( ^ ) X] /(^m,2/n)(sin27rA:x m sin27r/y n - sin 2yr/x m sin 2tt/c^ 



f m,n=0 1 
\ m>n J 



3.2.3. Example of antisymmetric interpolation. 

Consider the following multiple of the Gaussian distribution 



/(z,2/) = e" 



(51) 



(52) 

where (x',?/) = (0.707,0.293) and a = 0.079. The function /, restricted to the domain F(Sf ff ), is 
depicted in Figure |3j 





Figure 3. The function /, given by (52), plotted over the domain F(Sf). 



We sample the function / on the grids L 1 , L 1 and L 1 and calculate the antisymmetric 



interpolation functions ^ and tp± 2 fro m the relations (34) and (37), (38). These interpolating 
functions, together with the interpolating grids, are depicted in Figure [4j Interpolation errors are 
summarized in Table [2 

3.3. Symmetric interpolation. 

For interpolation with symmetric exponential functions, we consider the 'closed' triangle ^Ja+i] 
inside the square K^ a ^ : 



K [a,*+i] = 2/) ^ [ a 5 a + !] X [ a 5 a + !] \ x >y}- 



(53) 



For a given function / : ^J a+1 j —> C and a set of points L^ bN1 C ^J a +ip we define an symmetric 
interpolating function 

M 

^n( x iV)= Yl c ki E ( k ,i)( x ^)^ x,yeR (54) 

(k,l=-M\ 
I k>l J 



satisfying 



^(^m, J/n) J/n), TO > ft, TO, n = . . . TV — 1. 



(55) 
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AT = 4 



N = 7 



N = 12 






Figure 4. The antisymmetric interpolating trigonometric functions ^ 4 , tp 7 and 
of /, given by (52). The points of the interpolation grids L~ k , L~ k and L~ ± 

are depicted as small black dots. 



For TV = 2M+1, the antisymmetric interpolating polynomial has N(N+l)/2 unknown coefficients 



Cku to which correspond N{N +l)/2 constraints (55). For TV = 2M, we assume further 2M+1 conditions 

(56) 



t nCm,v l = -M...M 



and we have (2M + 1)(2M + 2)/2 - (2M + 1) = iV(iV + l)/2 unknown coefficients c\ { corresponding 
to N(N + l)/2 constraints (55). For the function / : ^^a+i] ~^ we define the symmetric extension 
Sf : i^[a 5 a+i] ^ C by the formula 



Sf(x,y) 



f(x,y) x>y 
f(y,x) x<y. 



(57) 



Using symmetric extension, one can prove the following proposition, similarly to Proposition 3.1 



Proposition 3.2. There exists a unique symmetric interpolating polynomial (54) satisfying (55). The 
coefficients are given for N = 2M + 1 by 



N-l 



°kl~Q jy2 5^ Grrmf( x miyn)Efai)(x m ,y n ) 
kl f m,n=0 | 
\ m>n J 



(58) 



and for iV = 2M assuming ( 56 ) by 



gfc.Mfl,M V G' 1 fix v )E + (x v ) 

kl r m,n=0 1 

I m>n J 



AT-1 



-kl 



(59) 



3.3.1. Calculation of the coefficients c 

Instead of the direct calculation of coefficients c^, the symmetric discrete Fourier transform (|22|) can 
be used, and the resulting coefficients transformed to c^'s. By direct comparison of (58), (59) to 



(22), we obtain for N = 2M + 1 



-kl ~ ^kl 



/?+, k>l,k,l = 0...M 



"t,-i = TN PN-i,k> fc = 0...M,Z = l...M 



(60) 



-l = T 2 N P 



N—k,N- 



.,, k < I, k,l = 1 . . . M 
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and for N = 2M 

4i = 9k,M9l,Mpki, 
c t-l = 9kM T N^N-i^ 
C -k,-l T N^N-k,N-V 



k > /, k,l = 0...M 

fc = 0. ..M, Z = 1...M- 1 (61) 
fc < Z, fc,Z = 1...M- 1 



The formula (61) determines (2M)(2M + l)/2 coefficients; the rest of the 2M + 1 coefficients are 
determined via the relations (l36l). 



3.3.2. Trigonometric form of 



The symmetric interpolating polynomial (54) can be brought to its trigonometric form: 



M 



ip~^(x,y) = hkhiG kJ l [A^ l (cos27vkx cos 27vly + cos2tvIx cos27vky) 

(k,l=0\ 
I k>l J 

+£>^ (sin 27r/cx cos 2nly + cos 27rZx sin 2nky) + 
+C^(cos 27r/cx sin 2nly + sin 27r/x cos 2ixky) 
(sin 27r/cx sin 2nly + sin 27r/x sin 2nky)] 



(62) 



where 



fc/ k k kj I kj I 

B kl =i ( c fc,/ ~ C -kJ + C fc,-/ ~~ C -fc,-/) 



(63) 



ki 



k k kj I kj I* 

k k k^ / k , / 



Substituting equations (37), (38) into (49) we obtain the following explicit formulas for N — 2M + 1 



^ki 



C ki 



N-l 



GkiN 2 



^ G r J n f(x m , y n )(cos 2nkx m cos 2nly n + cos 2nlx m cos 2nky n 

f m,n=0 "1 

4 N ~ 1 

B kl = r N 2 ^ G^/(x m ,y n )(sin27rfcx m cos27r/y n + cos27r/x m sin27rA;y n ) 

^ f m,n=0 1 
1 m>n J 



JV-1 



^ ^ ^ G m ^/(x m ,y n )(cos27r/cx m sin27r/y n + sin 2yr/x m cos 2tt/c^ 

^ f ra,n=0 1 
1 m>n J 



JV-1 



^ifei = r ^ G m ^/(x m ,y n )(sin27rfcx m sin27r/y ri + sin 2yr/x m sin 27rfcy n ) 

^ r ra,n=0 1 
1 m>n J 



(64) 
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and for N = 2M 



A kl = 4 ^ M jv2 M G mnf( x m, y n )(cos 27rfcx m cos 27r/y n + cos 2yr/x m COS 27T% n ) 

\ m>n J 



AT-1 



5 



^9k,M9l,M 



N-l 



ki - r AT2 yZ Gj n f(x m ,y n )(sm27rkx m cos27rly n + cos27rlx m sm27rky n ) 

\ m>n J 

7 2 ^ G^/(x m ,^ n )(cos27r/cx m sin27r/y n + sin27r/x m cos27rA:y n ) 

^ f m,n=0 1 

\ m>n J 



^9k,M9l,M 

GmN 2 



N-l 



^2 G m 1 n f(x m ,y n )(sm27rkx m sin 2nly n + sin 27vlx m sin 2nky ri 

( m,n=0 "1 
\ m>n J 



(65) 



3.3.3. Example of symmetric interpolation. 

We sample the function /, given by (52), on the grids , Lj" x and , and calculate 

U, 9,4,1 (J, 7>, 7,1 U,7J,1^,1 



the symmetric interpolation functions ^ , ^ and fro m the relations (54) and (|58| ), (59). These 
interpolating functions, together with the interpolating grids, are depicted in Figure |5J Interpolations 
errors are summarized in Table [TJ 






N = 4 



N = 7 



N = 12 






Figure 5. The symmetric interpolating trigonometric functions tp^ and ^)\ 2 of /, 
given by (52). The points of the interpolation grids L+i 4 ^ ^1 7 x an d Z^J~i 12 1 are 

depicted as small black dots. 



4. ( Antisymmetric cosine transforms 
4.1. Antisymmetric cosine functions. 

4.1.1. Definitions, symmetries and general properties. 

Two-dimensional antisymmetric cosine functions cos^ ^ : R 2 — > C have for A,/i G K the following 
explicit form 

cos" M) (:r, y) = 1^(1^ = cos(^Ax) cos(tt W ) - cos^/xx) cos(^Ay) 



COs(7T/J J x) COs(7T/iy) 



Note that we chose this definition according to Section X of [7J . Instead of the factor 2tt we use 'half 
argument 7i\We observe that cos^ A ^(x,y) = and cos^^(x,x) = 0. From the explicit formula, we 
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immediately obtain antisymmetry of cos^ ^(x, y) with respect to the permutation of variables (x,y) 

COS (A,m)^' X ) = - COS (A,/x)^'^- ^ 66 ) 

and, moreover, with respect to the permutation of (A, /i) 

COS (A,/i)^'^ = ~ COS (m,A)( X '^)' ( 67 ) 

Therefore, we consider only such cos^^ with strictly dominant (A,/x), A > /i. The functions cos^.^ 
with fc,lEZ have symmetries related to the periodicity of cosine function 

cos^^ (x + 2r,y + 2s) = cos^ z) (x,y), r, s E Z 
We also have invariance under the change of sign of (x, y) 

cos ( k ,i)(- x >y) = cos (fc,o^' = cos (m)(~ x ' = cos (k,i)( x >y) 

and invariance under the change of sign of (fc, /) 

cos ("-fc,o = cos (fc -/) = cos ("-/c -o ( x ' y) = cos (kj) ( x ' y) • ( 70 ) 



(68) 
(69) 



The relations (|66|) - (|70|) imply that it is sufficient to consider the functions cos^^, k,l E Z-°, fc > / 



on the fundamental domain F(Sf). 

4.1.2. Continuous orthogonality. 

The functions cos^, ^ are mutually orthogonal on the fundamental domain F(S^), i.e., 



/ cos (k n (x,y)cos fife/ „Jx,y)dxdy. = ^5 kk t5 w , k,l,k\l 

JF(Sf) V ' ,v V '' '' 4 



E Z^ u , k > Lk > I 



Also every function / : R 2 — ► C that is antisymmetric f(x, y) = — /(y, and periodic f(x + 2r, y + 2s) = 
/(#, y), r, 5 E Z and has continuous derivatives can be expanded in the antisymmetric cosine functions 
cos7 



f(x,y) = ^2 c fc/cos (/c/) (x,y), c H = 4:J^^f(x,y)cos^(x,y)dxdy 



k>i 



(71) 



The graphs of the lowest antisymmetric cosine functions cos^ ^, k, I E {0, . . . , 3}, k > I are plotted 
in Figure [6| 




^(2,1) 




°(3,1) 




(3,2) 



Figure 6. The contour plots of the antisymmetric cosine functions cos^^, k,l E 



{0,...,3}, k > 1. 
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4.1.3. Solutions of the Laplace equation. 

The functions cos - are solutions of the Laplace equation 



( d 2 d 2 



(72) 



and moreover of the equation 



d 2 d 2 



cos,, n (x,y) = 7T 4 fc 2 / 2 cos,, n (x,y). 



dx 2 dy 2 OMr > y) ~ ^{k,iy 
The functions cos - satisfy the condition cos^, ^(x, y) — on the boundary x — y and the condition 

d 

— cos m (x,y) = 0, 
where n is the normal to the boundary x = 1 or y = 0. 

4.2. Antisymmetric discrete cosine transforms. 

Four types of discrete antisymmetric cosine transforms can be derived from the antisymmetric trigono- 
metric transform. In order to derive these transforms, we define the following two functional operators. 
To a function / : i^[o,L] we assign a function E^f : i^[_L,L] ~~ > C defined by 



E L f(x,y) :-- 



'f(x,y) x>0,y>0 

f(-x,y) x<0,y>0 

f(x,-y) x>0,y<0 

J(-x,-y) x<0,y<0, 



(73) 



and to a function / : if[o,i] ^ C, we assign a function Rf : K [ 0j 2] — > C defined by 

f/(x,y) 0<^< 1,0 < 2/ < 1 

-f(2-x,y) Kx<2,0<y<l 

-f(x,2-y) 0<x<l,Ky<2 

U(2-x,2-y) Kx<2,Ky<2. 



Rf(x,y) 



(74) 



All four antisymmetric cosine transforms operate on a function / : ^ — > C. Applying the formula 
5ll) to the four functions 



(I) E x Af : K[_ M] C, where N = 2M, T = 2, 6 = 1 
(II) A/ : K[_ M] C, where TV = 2M, T = 2, 6 = 1/2 

(III) £ 2 iM / : K h252] -> C, where N = 4M, T = 4, 6 = 1 

(IV) E 2 RAf : ^[-2^2] -> C, where TV = 4M. T = 4, 6 = 1/2 

we obtain due to (anti)symmetry in (48) that = = = 0. Introducing the symbol d^^M for 
fe = 0,...,Mby 



dk. 



M 



_\\ iffc = 0,M 
I 1 otherwise 



(75) 



the four interpolating functions (48) corresponding to cases (I)-(IV) and their coefficients (51) can be 
brought to the following form. 

AMDCT- I. 

M ^ ^ M 

^m"( X > V) = ^2 4',7 COS (fc,0^' ^' C M = M 2 ^ d m , M d n ,Mf (%, J/n) cos (fc,0 2/n) 

rfc,/=o\ ^ r m ,n=0 ], 



where x m = ^, y n = ^. 



M ' 
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AMDCT- II. 



M-l A r 7 M-l 



/II,-/ \ II- - / \ II- ^d k M di,M r ( \ - ( \ 

(x,y)= 2^ C M cos (M)^'^' C M = 1^ f(xm,Vn)coB M (x m ,y n ) 

(k,l=0\ f m,n=0 1 



AMDCT- III. 



, m+ - 77+ - 

where x m y n 



M-l 

*I>m~( x >v) = X] c S'" cos ^ + i^ + i)( x ^)' 

rjfe,z=o\ 2 ' 2 

I fe>Z / 



where x r 

AMDCT- IV. 



4 M - x 



f m,n=0 1 
\ m>n J 



ILL oi — — 
M ' W ~~ M " 



M-l , M-l 

4 



^M (^?/)= ^ 4T' COS ( fc+ i 5/+ i)( X ^)' 4T' =]^2 5^ /(Sm,!/n)coS (jfe+ i ^ (ff m , J/ n ) 

r fe,z=o i 2 2 r m^^oi 2 2 

\ k>l I I m>n J 



77+ 1 



where x m y n 



4.2.1. Example of antisymmetric cosine interpolation. 

We sample the function / on the grids L~i , L~ ± and L~i and calculate the antisymmet- 

ric cosine interpolating function of the type AMDCT-II: V^ 1 ' - , V^ 1 '" and V^' - - These interpolating 
functions, together with the interpolating grids, are depicted in Figure [7| Interpolation errors are 
summarized in Table [U 




Figure 7. The antisymmetric interpolating cosine functions of the type AMDCT-II: 
V^ 1 ' - ? ^7 I,_ an d ^12 ~ of the function /, given by ( [52] ) . The points of the interpolation 
grids L~ 1 , L~i and L~ ± are depicted as small black dots. 

5 O* 5 5 5 O* 5 5 5 0*5 5 



4.3. Symmetric cosine functions. 
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4.3.1. Definitions, symmetries and general properties. 

Two-dimensional symmetric cosine functions cos^ ^ : R 2 — > C have for A, fi E K the following explicit 
form 

cos+ (x, y) = = cos(^ Ax) cos(^) + cos(^x) cos(^Ay) 



COs(7T/ix) COs(7T/iy) 



Note that we chose this definition according to Section X of [7]. Instead of the factor 2tv we use 'half 

argument tt. From the explicit formula we immediately obtain symmetry of cos^^(x,y) with respect 
to the permutation of variables (x, y) 

cos (a, M )(^' x ) = cos (a, M )( x ' v)- ( 76 ) 

and, moreover, with respect to the permutation of (A, /i) 

cos+ M) (x, y) = cos J x) (x, y). (77) 

Therefore, we consider only such cos^^ with dominant (A,/x), A > \i. The functions cos J, ^ with 
GZ have symmetries related to the periodicity of cosine function 

cos J.^ (x + 2r, y + 2s) = cos^ z) (x, y) , r, s E Z. (78) 
We also have invariance under the change of sign of variables (x, y) 

COS J,0^~ X '^ = cos 0M)( x ' ~^ = cos (W~ x ' = cos Cfe,o( x ' y ) ( 79 ) 
and under the change of sign of (fc, /) 

cos ("-M)^'^ = cos J-/)(^2/) = cos^ fc _ z) (x,y) = cos+ M) (x,y). (80) 



The relations (|76|) - (|80|) imply that it is sufficient to consider the functions cosj, ^, fc,i G Z-°, fc > / 



on the fundamental domain F(Sf). 

4.3.2. Continuous orthogonality. 

The functions cos J, ^ are mutually orthogonal on the fundamental domain F(S^), i.e., 



Gki 



cost n (x,y)cost, v Ax,y)dxdy. = —Skk'Sw, k,l,k\l' E Z-°, k>l,k'> l'. 

Also every function / : R 2 — > C that is symmetric /(x,y) = /(y, and periodic /(x + 2r, y + 2s) = 
/(x,y), r, s E Z and has continuous derivatives can be expanded in the symmetric cosine functions 
cost 



f(x,y) 



V c fc/ cos+ (x,y), Cki^iG^ 1 f(x,y)cost kl) (x,y)dxdy. (81) 



Mez>°l - (^ 2 aff ) 



The graphs of the lowest symmetric cosine functions cos J, ^, fc, Z E {0, 1, 2}, fc > / are plotted in Figure|8j 

4.3.3. Solutions of the Laplace equation. 

The functions cos + are solutions of the Laplace equation 

and of the equation 

d^d^ cos (fe,0^'f) = vr 4 fc^ 2 cos+ M) (x,y). 
The functions cos + satisfy the condition 

^cos+ M) (x,y) = 0, 

where n is the normal to the boundary of F^SqP), i.e. are the solutions of the Neumann boundary value 
problem for F(S^). 
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(0,0) 













m 


COS^ 


1,1) 




°(2,1) 




(2,2) 



Figure 8. The contour plots of the symmetric cosine functions cosj,^, k,l G 



{0,1,2}, k > I 



4.4. Symmetric discrete cosine transforms. 

Four types of discrete symmetric cosine transforms can be derived from the symmetric trigonometric 
transform. All four symmetric cosine transforms operate on a function / : — > C. Applying the 
formula (65) to the four functions 



(I) EiSf : X[_ M] C, where N = 2M, T = 2, b = 1 
(II) ExSf : K[_ M] C, where TV = 2M, T = 2, 6 = 1/2 

(III) E 2 RSf : #[12,2] -> C, where N = 4M, T = 4, 6 = 1 

(IV) £ 2 #S/ : K[_ 2 ^ 2] -> C, where N = 4M, T = 4, 6 = 1/2 



we obtain due to symmetry in (62) that = = 0. The four interpolating functions (62) 



corresponding to cases (I)-(IV) and their coefficients (65) can be brought to the following form. 
SMDCT-I. 



M 



y)= J2 ° 1 m cos tk,i)( x > y)' 

r k,i=o \ 
I k>l J 



i,+ _ 4dfe,M^,M r^- 1 ^ a *( \ + ( 

c kl ~ M 2 G J ^ r mn a m 1 M^n,M J {%m, Vn) cos (& ? j) l x m 5 ; 



7U 



f m,n=0 "1 
I m>n J 

where x — m — n 



SMDCT-II. 



m M' 



M ' 



M-l 



M-l 



^M + (^> ?/) = 51 C £ +COS J,o( X ' C £ + = M ^tn M G mlif ( X ™, Vn) COsJ^ (x m , y n 

r 7„ l n >, kl r _ n >. 



I /c>£ J 

where x 



f m,n=0 1 
I m>n J 



ra+7 



SMDCT-III. 



m — M ' 



J/n 



M 



M-l 



rjfe,z=o\ 2 ' 2 

I fc>Z j 



M-l 



III, 



M 2 G 



ki 



51 dm,M d n,M G mnf ( x m, Vn) COS J + i >z+ i) ( 



f m,n=0 1 
\ m>n J 

where x m = ^, y n = ^. 
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N 



J F (Sf) ~ f[ 



F(S: 



/II,- 



f 



F(Sf) 



,II,H 



f 



97987 



10" 



97336 • 10 



94170 



io- 



89002 



10" 



86234 



10" 



86224 



10" 



77865 



10" 



77839 



10" 



21116 



10" 



21447 



10" 



35708 



10" 



10^ 
10^ 



35636 



10" 



10^ 
10^ 



9841 



10 



9812 



10" 



14023 



13915 



1949 



10 



1978 



10 



2570 



10=^ 



2570 



10=^ 



1000 



10 



1001 



10 



1309 



1310 



10 



503 



10 



504 



10 



600 



601 



11 



63 



10" 



63 



10 



86 



10" 



86 



10" 



12 



10" 



10" 



11 



io- 



ii 



10" 



Table 1. Comparison of errors of interpolations %p N , ip^, ip 1 ^' 
/, given by (52), is sampled on grids L^ x , N = 4, ... ,12. 



and ip 1 ^' 



The function 



SMDCT-IV. 



M-l 



/ j ^kl 
(k,l=0\ 
X k>l J 



M-l 



„IV,H 
-Jfci 



M 2 G kl 



^1 Gmnf ( x mi Vn) cos J, + i { x rm V\ 

( m,n=0 1 
\ m>n J 



wnere x m — ^ , y n — ~~m~ • 
4.4.1. Two Examples. 

(1) We sample the function / on the grids 1 , 1 and 1 and calculate the symmetric co- 

2 ' ' ' 2 ' ' ' 2 ' ' 

sine interpolating function of type SMDCT-II: f/^ 1 ^, V*?^ an d ^i2 + - These interpolating functions, 
together with the interpolating grids, are depicted in Figure |9j Interpolation errors are summarized 
in Table [Q 






M = 4 



M = 7 



M = 12 






Figure 9. The symmetric interpolating cosine functions of type SMDCT-II: ^"' + , iff*' 



of /, given by (52). The points of the interpolation grids L 



and ^2 

. are depicted as small black dots. 

U, g ,1^,1 



o,Ui' 



L 0,i,7,l aM 
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(2) Symmetric cosine interpolating functions ip 1 ^ and ^^' + offer an additional advantage due to the 
absence of the so called Gibbs phenomenon alongside the borders of F(Sf ff ). The Gibbs phenome- 
non for one-dimensional Fourier expansions, which describes specific behavior ('ringing') of Fourier 
expansions at the points of discontinuities of / was extensively studied, see e.g. the detailed re- 
view [5j . The Gibbs phenomenon in dimension two was also investigated in \17\ H] . The analog of 
the Gibbs phenomenon appears in Fourier interpolating functions [3J. Since the (anti)symmetric 
and periodic extensions of continuous functions on F(*S| ff ) may generate discontinuities of the re- 
sulting functions on the borders of F(*S| ff ), the Gibbs phenomenon is in general to be expected for 
the interpolating functions 

• ifj^: on the borders y = 0, x = 1 and x = y 

• tp^: on the borders y = 0, x = 1 



and ^j l b : on the border x 



Since for the continuous function f on F(Sf ff ), the periodic extension of E\Sf (see (57), (73)) is 
continuous, the Gibbs phenomenon for the interpolations and of f r does not occur, 
fact is illustrated in Figure [TO} 



his 





M = 12 





M = 20 




Figure 10. The symmetric interpolating cosine functions of type SMDCT-II: ^ 
^2o + of the function f(x, y) = ^ cos 40^ = 0.707) 2 + (y - 0.293) 2 . 



5. Concluding remarks 

When functions of only two variables are considered, the structure of our special functions is quire 
transparent. For this reason, we avoided any reference to the underlying symmetric group S3 of per- 
mutations of three elements [T6] . 

Section 4 contains a description of symmetric and antisymmetric 2D transforms built on cosine 
functions. An analogous presentation of the properties of the transforms built on sine functions would 
be of interest. 

None of the six families of special functions considered in the paper correspond to the well known 
cosine or sine transforms [T51 fT3] , where one dimensional transforms are used simultaneously in two mu- 
tually orthogonal directions. Nevertheless it is justifiable to claim that here we have 2D generalizations 
of common trigonometric functions of one dimension. It is curious to see that both sine and cosine each 
have a symmetric and antisymmetric generalization. 

The conventional cosine transform used in 2D comes in four types |18| , traditionally labeled I, II, III, 
and IV. They differ by the shift of the data sampling points with respect to the original lattice points. 
The shifting is by a half distance separating the lattice points in either of the two space directions or 
in both. 

A comparison of theoretical properties of the 'cosine transforms of S3' in this paper and those of the 
standard S2 x S2 would be of interest and has yet to be made. The arguments in favor of the S3 version 
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of the transforms quote the ease of generalization to any dimension [6J and the possibility to work with 
data on lattices of other symmetries |16| . In general the greater the symmetry group underlying the 
formalism, the more economies one may expect in some applications. On the down side, there is the 
need to work with functions situated in domains that are of triangular shape even in square lattices. 

Practical computational aspects of the formalism presented in this paper need further investigation. 
The identification of problems in which it can be most advantageous, a comparison of computation 
speeds, and modification to the fast Fourier transform, etc. should be addressed. 

Visual inspection of the interpolations of the Gaussian model functions in this paper lead to the 
qualitative conclusion that the interpolation error is rather small once the minimal distances of the 
lattice grid become smaller than the dispersion of the Gaussian model function. 

The orbit functions of the paper have other useful properties that were not exposed here, two of which 
are: (i) Their products are decomposable into their finite sums, (ii) They can be naturally rewritten in 
terms of variables referring to non-orthogonal bases related to the simple roots of the simple Lie group 
SU(3), see [16]. 
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